Thermal conductivity of pure and Mg-doped CuGeC>3 in the incommensurate phase 
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The thermal conductivity k, of Cui-^Mg^GeOa is measured in magnetic fields up to 16 T. At 
the transition field H c to the high-field incommensurate (I) phase, k abruptly decreases. While 
k of the pure CuGeC>3 is enhanced with the application of higher magnetic fields, an anomalous 
plateau feature shows up in the k(H) profile in the I phase of the Mg-doped samples which includes 
antiferromagnetic ordering (I-AF phase). With the help of specific heat data, taken supplementally 
for the identical samples, the above features are well understood that phonons are strongly scattered 
by spin solitons and that n in the I phase is governed by their spatial distribution. In particular, the 
plateau feature in the doped samples suggests "freezing" of the solitons, in which antiferromagnetism 
must be playing an essential role because this feature is absent at temperatures above Tm- 

PACS numbers: 66.70.+f, 75.30. Kz, 75.50.Ee 



I. INTRODUCTION 

Recently, a new class of commensurate- 
incommensurate (C-I) transition, magnetic-field induced 
C-I transition, which had been theoretically predicted,El 
was experimentally established in an one-dimensional an- 
tiferromagnetic (AF) Heisenberg spin system coupled to 
lattice degrees of freedom, the spin-Peierls (SP) system; 
the magnetic n tcansition is observed in magnetization 
measurements&tj and the incommensurate lattice mod-it. 
lation is visibly shown in x-ray diffraction experiments qEI 
The SP system shows a transition at Tsp from the uni- 
form (U) spin- 1/2 quantum liquid phase to the dimerized 
(D) phase, where spin-singlet ground state is focmed at 
the expense of the lattice deformation energy.H With 
the application of magnetic field, a new magnetic phase 
appears above a threshold field H c because the nonmag- 
netic spin-singlet state becomes energetically unfavorable 
in high magnetic fields. The magnetic moment is partly 
recovered in the chain, and "spin solitons" are formed 
in this magnetic phase. The "spin soliton" contains 
both lattice deformation and moment (carrying spin 
1/2), whose spatial distribution is governed by the sine- 
Gordon equation. Usually, the periodicity of the solitons 
are determined simply by the applied magnetic field, and 
as a result, the spin and lattice modulation becomes in- 
commensurate. Therefore, this high-field phase is called 
incommensurate (I) phase. After the existence of the in- 
commensurate modulation is established, the subsequent 
interest on this I phase is how the spin and lattice mod- 



ulation evolves under higher magnetic fields (well above 
H c ), where the soliton distance becomes comparable to 
the soliton width. LTQ 

The discovery of the inorganic SP compound CuGeC>3 
(Ref. |l0| ) opened the way to impurity substitiution and it 
turned out that just a small amount of the impurity mod-. 
ifi.es the phase diagram to a more complicated oneJiMHl 
with decreasing temperature, after the appearance of D 
phase at Tgp, the system shows another transition to 
a peculiar ordered phase, where antiferromagnetic (AF) 
long-range ordei^and the SP order coexist [dimerized AF 
(D-AF) phase] £3 The staggered moment and the lattice 
deformation are spatially modulated in D-AF phase; the 
amplitude is larger, .near the impurity sites and smaller in 
between themoo In contrast to such well-established 
pictures of the ordered phase in low fields (D and D-AF 
phase), the spin distribution in high-field phases of the 
impurity-doped systems is still under question. Also, the 
phase boundaries for the high-field phases are not yet 
investigated for the system with its Cu site substituted 
by impurities (Cu-0 chains are directly disordered), al- 
though there exist several reports on these phase bound- 
aries for the system whose Ge site is substituted. Re- 
cently, we reported that k of Cu^rMg^GeOs is sensi- 
tive to the Mg-doping (disorder) JlaEZl and demonstrated 
that such transport measurement is useful in examining 
the disorder effects. In this work, we have measured the 
thermal conductivity of Cui-^Mg^GeOa single crystals 
with the application of magnetic field up to 16 T. The 
specific heat is also measured supplementally for the iden- 
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tical samples. First, k of the pure CuGeC>3 is measured in 
order to identify the dominant scattering mechanism in 
the I phase. It turned out that the dominant heat carri- 
ers (phonons) are strongly scattered with the appearance 
of the spin solitons in the I phase, indicating that k is a 
good probe of the soliton distribution. Our main result 
is an anomalous plateau feature found in the k(H) pro- 
file in the antiferromagnetically ordered I (I-AF) phase of 
the Mg-doped samples, suggesting that the spin solitons 
are frozen in the FAF phase. 



II. EXPERIMENTAL 



The single crystals of Cui_ x Mg x Ge03 were grown with 
a floating-zone method. We have prepared both pure 
CuGe03 (the SP transition temperature Tsp=14.5 K) 
and two Mg-doped samples x = 0.016 and x = 0.0216 
with T S p's of 11.5 K and 9 K, respectively. (The T SP 's are 
determined by dc magnetization measurements) . The Mg 
concentration is carefully determined by inductively cotu 
pled plasma-atomic emission spectroscopy (ICP-AES).t3 
The Mg-doped samples shows the .-Keel transition at 
T N = 2.5 K and 3 K, respectivelyJi3'E3 Masuda et al. 
carefully examined x dependence of the TV for the 
same series of crystals and found that TV jumps at the 
impurity-driven transition from the D-AE. phase to the 
uniform AF phase without dimerization£3'E£l Note that 
the D-AF phase shows up below in the sample used 
for our present study because the Mg concentration of the 
both samples is smaller than the critical concentration x c 
(= 0.023 - 0.027) Ma 

The thermal conductivity is measured using a "one 
heater, two thermometers" method. The direction of the 
heat current and the magnetic field is parallel to the o 
axis, i.e., the chain direction. The detailed experimen- 
tal technique is described elsewhere .E3E3 Since the total 
thermal conductivity is expressed as k = ^CiDi, where 
Ci and Di are specific heat and diffusivity of each heat 
carrier denoted by i (the diffusivity Di is proportional to 
the mean free path of the heat carrier), the specific heat 
measurement is helpful in the analysis of the n data. We 
have measured the specific heat in magetic fields up to 16 
T by the usual relaxation method for the same samples. 
The addenda consists of a sapphire substrate, a cernox 
temperature sensor (calibrated in magnetic fields) and 1- 
kf2 micro heater. The connection from the addenda to 
the heat sink is made by Sb-Au thin wires to provide a 
moderate heat-relaxation rate, from which the heat ca- 
pacity is calculated. 



III. RESULTS 



A. Pure CuGe0 3 
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FIG. 1. (a) Magnetic field dependence of the thermal con- 
ductivity of CuGeC>3 single crystal at 4.2 K. The magnetic 
field is applied in the c-axis direction. First order transi- 
tion from the spin-Peierls to the incommensurate phase takes 
place at H c . Inset: temperature dependence of the thermal 
conductivity of CuGeC>3 in magnetic fields. Tsp and Tj are 
the transition temperatures below which spin-Peierls phase 
and incommensurate phase appear, respectively, (b) C'/T in 
magnetic fields is plotted as a function of T 2 . The dotted line 
for the 0-T curve indicates the T 3 dependence of the phonon 
contribution C p h- The other dotted line for the C/T curves in 
the I phase shows that the additional magnetic contribution 
C m also depends on temperature as T 3 . 

The inset of Fig. 1(a) shows the temperature depen- 
dence of k in fields up to 16 T. The transition tempera- 
tures from the U phase to the D phase (below 9 T) and 
to the I phase (above 14 T) are indicated as Tsp and Tj, 
respectively. While both phonons and spin,fixcitations 
carry a large amount of heat in the U phase,E3Ej phonons 
become the dominant heat carrier at temperatures well 
below Tsp in the D phase because of the spin gap. These 
phonons are mostly scattered by defects at the lowest 
temperatures, where spin excitations are almost absent, 
and field dependence is small there; the 0-T and 9-T 
curves merge below ~ 3 K. When thermally excited spin 
excitations increase with temperature and the scattering 
rate of the phonons are enhanced, via spin-phonon cou- 
pling, k starts to decrease with temperature, producing 
a pronounced peak around ~ 5 KE3 Since this peak in 
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k is a manifestation of the spin gap, the-.peak height is 
rapidly suppressed with magnetic fields J13 

The field dependence of k at 4.2 K is presented in the 
main panel of Fig. 1 (a) to see the suppression of the peak 
under magnetic fields. A sudden drop in n, accompanied 
by hysteresis, shows up at the first order transition from 
the D phase to the I phase (H c ~ 12.5 T). When the mag- 
netic field increases above H c , n starts to increase at ~ 14 
T. The specific heat C is measured for the same sample 
in order to elucidate the origin of the field dependence 
of k {Ki = CiDi). The 14-T and 16-T data are plotted 
together with zero-field data in Fig. 1(b). In contrast to 
k, C in the I phase is larger than that in the D phase 
below Tj ~ 10 K, indicating that the jump in k at H c is 
caused by a sudden change in diffusivity. It is also to be 
noted that the values of C at 14 T and at 16 T do not 
significantly differ below Tj, while n shows an apparent 
increase above 14 T, suggesting the enhancement in the 
diffusivity Di with the field application above ~ 14 T. 



B. Mg-doped CuGe0 3 
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FIG. 2. (a) Magnetic field (along the CuO chain) depen- 
dence of the thermal conductivity of Cui-zMg^GeOa with 
x = 0.016 at 5 K. (b) Temperature dependence of k in mag- 
netic fields up to 16 T for the x = 0.016 sample. Tjv indicates 
the Neel transition in the I phase. 

Figure 2(a) shows the k(H) profile of the x = 0.016 
sample at 5 K, which is above the zero-field (~ 2 K), 
determined from dc magnetization. Similarly to the pure 
sample [Fig. 1(a)], k is rapidly suppressed with magnetic 



field when the field approaches H c ~ 11.5 T, and then 
increases in the field range above H c . Hysteresis is ob- 
served in the vicinity of H c also for this sample, although 
the transition becomes broader. 



The temperature dependence of k is shown in Fig. 
2(b). The transition temperatures Tsp an d Tj are in- 
dicated with arrows (together with T/y). Note that the 
14-T and 16-T curves merge below Tjv, suggesting field 
independent k(H) above H c in this temperature range. 
Since this feature is expected to be more pronounced for 
a sample with higher Tjv, we measured the field depen- 
dence of the x = 0.0216 sample. Figure 3(a) presents 
n(H) profiles at the temperatures of 2.5 K and 5 K. An 
apparent plateau is observed above H r ~ 11.5 T in the 
2.5-K curve, while the 5-K curve shows a slight upturn 
above H c , similarly to the field profile in the I phase of 
the pure sample. Figure 3(b) shows the temperature de- 
pendence of k in three different fields above H c and that 
without magnetic field. One can see that the 12-, 14- 
and 16-T curves merge below ~ 4 K, which is close to 
zero-field T N (~3K). 
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FIG. 3. (a) Magnetic field dependence of the thermal con- 
ductivity of Cui-zMg^GeOa with x = 0.0216 at 2.5 K. k(H) 
becomes independent of magnetic fields in the I phase (above 
H c ). (b) k in magnetic fields up to 16 T is plotted against T 2 
for the x — 0.0216 sample. Below Tjv, all the curves for the 
I phase merges to a single line. The dotted line shows that k, 
is mostly proportional to T 2 below Tjv. 

Figure 4 shows the temperature dependence of the spe- 
cific heat of the same x = 0.0216 sample. The 16- T curve 
shows a A-shaped anomaly, indicating that Tjv under this 
magnetic field is around 4 K, just below which the plateau 
feature appears in the n(H) profile. Note that Tjv at 16 T 
is apparently higher than the zero-field T^r, because the 
spin fluctuation is reduced by field applicationE3€il It is 
also notable that a broad anomaly is observed above Tn 
up to ~ 8.5 K, suggesting the existence of another phase 
at higher temperatures. 
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FIG. 4. Specific heat of the x = 0.0216 sample. The I 
phase and the I-AF phase appear below Tj and Tn, respec- 
tively. The solid line is a guide to determine Ti. 



IV. DISCUSSION 
A. Soliton scattering in I phase 

First, let us consider the k(H) profile at 4.2 K shown in 
Fig. 1 (a) . The dominant heat carriers are phonons in the 
pure CuGeC>3 well below Tsp, where the spin excitations 
are hardly excited because of the spin gap. Therefore, 
k(H) shown in Fig. 1(a) mostly reflects the field depen- 
dence of K p h at least below H c . When the population 
of the spin excitations gradually increases with field, the 
scattering by the spin excitations increases and, as the 
result, K p h slowly decreases. 

While the major part of the specific heat is due to 
phonons in the D phase, the additional component, i.e., 
magnetic specific heat C m , appears in the I phase, as 
shown in Fig_l(b). This result reproduces the report by 
Lorenz et a/.,E3 and is consistent with the theory claim- 
ing that_C m is due to magnetoelastic excitations called 
phasons a Since both C m and C p h depends on temper- 
ature as oc T 3 , the total C is cx T 3 . The thermal con- 
ductivity is given by k — ^CjDj (i = ph, to), where 



the additional magnetic specific heat C m reaches twice 
as large as the phonon contribution C p h [see Fig. 1(b)]. 
Nevertheless, k suddenly decreases at the transition from 
the D phase to the I phase, meaning that the decrease 
in K p h is dominant at the transition, though the newly 
appearing phasons could also carry some heat. Since lit- 
tle change in C p h is expected at H c , it must be -D p h that 
is responsible for the abrupt change in k(H), indicating 
the sudden appearance of new scatterers of phonons in 
the I phase. It is most likely that phonons are strongly 
scattered by domain walls accompanied with the spin 
solitons, which would be equivalent to the scattering by 
the phasons. Note that the ultra-sonic results, in which 
sudden change in elastic constants was shown at H c , axe 
understood in terms of the phonon-soliton interaction.cil 

The k(H) profile above H c also is not simply explained 
by C in the I phase; K increases above ~ 14 T, while the 
difference in C m at 14 T and that at 16 T is negligi- 
ble within the accuracy of our measurement (~ 10% of 
C m )j23 as shown in Fig. 1(b). Since the enhancement 
in k from 14 T to 16 T is more than four times larger 
than this maximum error, the diffusivity should be re- 
sponsible for the magnetic-field dependence of k also in 
the field range above H c . 

Assuming that D p h(H) is governed by the scattering 
by the domain walls in the I phase, as is discussed above, 
Dph(H) should be modified when the soliton distribu- 
tion changes with field. In a sufficiently high magnetic 
field where the soliton distance d becomes comparable 
to the soliton width £, the adjacent solitons would over- 
lap with each other and the lattice and spin modulation 
might become sinusoidal wave, as is suggested by a recent 
NMR measurementQEl and a magnetostriction result .□ It 
is expected that the cross section of the domain- wall scat- 
tering is smaller for the sinusoidal-wave modulation than 
for the independent solitons, because the amplitude of 
the lattice modulation is smaller for the sinusoidal-wave 
distribution. Such field dependence of the soliton distri- 
bution well explains the enhancement in D p h(H) with 
field above H c , and is consistent with the k(H) profile in 
the I phase. The ratio of d/£ is - 2.3 at ~ 14 T (Ref. |) 
where the overlapping becomes evident in k(H). 



B. Incommensurate phases in the Mg-doped samples 



Since k and C data demonstrate anomalies at the tran- 
sition temperatures and at the critical field, we can dis- 
cuss the phase boundary for the high-field phases of the 
Mg-doped CuGeC>3. So far, contradictory results are 
reported on the high-field incommensurate phases for 
CuGe03 crystals whose Ge site is substituted by Si. Two 
transitions are observed in the specific heat measurement; 
one is from the U phase to the incommensurate phase 
without AF ordering (I phase) and the otherjs a Neel 
transition from this I phase to the I-AF phased On the 
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other hand, the result of a ultrasound-velocity measure- 
ment suggests a direct transition from the U phase to the 
I-AF phaseEa 

Let us discuss this problem on the basis of our k and C 
data for the Mg-doped samples, where Cu-0 chain is di- 
rectly modulated. The Neel transition to the I-AF phase 
is apparent also for the Mg-doped samples; the specific 
heat shows a clear A-shaped peak (Fig. 4), and the field 
dependence of k disappears as shown in Fig. 2(b) and 
Fig. 3(b). In addition to the features at the Neel order- 
ing, an enhancement is observed below Tj ~ 9 K in k 
for the x = 0.016 sample [Fig. 2(b)]. The specific heat 
data for the x = 0.0216 sample in 16 T also shows an 
enhancement below Tj ~ 8.5 K (Fig. 4). The results 
indicate the existence of the, I phase without AF order 
for the Mg-doped samples. EZJ However, the anomaly in 
the C data is not so obvious at Tj as at Tjv. This is 
probably because of the existence of short-range ordered 
(SRO) state in the SP transition in the impurity-doped 
CuGeC>3. Noting that long-range SP order is not estab- 
lished down to ~ 6.5 K at T for the x = 0.0216 sample 
(Rcf. 18 |2S|) , and assuming the SRO state also in mag- 
netic fields above H c , the "transition temperature Tj" 
will not be well-defined. 




N 

FIG. 5. The schematic H-T phase diagram of the 
x = 0.0216 sample, based on the results of the k and C 
measurements. Open and closed circles indicate the transi- 
tion temperatures determined from the k and C data. The 
hatched area indicates the short-range ordered regime. 

The field induced transition from D to I and from D- 
AF to I-AF phases at H c is also evident in the n(H) 
profile; K for the x = 0.016 sample and the two curves 
for the x = 0.0216 sample (at 2.5 K and 5 K) rapidly 
decrease (with hysteresis when x = 0.016) with the ap- 
plication of magnetic fields below H c ~ 11.5 T, while 
the decrease abruptly disappears above H c [Fig. 2(a) 
and Fig. 3(a)], similarly to the D-I transition observed 
for the pure CuGeC>3. Note that the transition at 5 K 
(above T N ~ 4 K) and that at 2.5 K (below T N ) take 



place at almost the same field around H c ~ 11.5 T, i.e., 
the value of H c is irrelevant to whether the temperature 
is above or below Tjv. The above results suggest that 
both transitions, from D to I and from D-AF to I-AF, 
are "commensurate-incommensurate" transition, which 
is characterized by the appearance of the spin solitons. 
The suggested phase diagram is presented in Fig. 5, fol- 
lowing the discussion in this subsection. 



C. Freezing of the solitons in I-AF phase 

The most striking observation in the k(H) profile for 
the Mg-doped samples is the plateau feature which ap- 
pears in the I-AF phase [see Fig. 3(a)]. Since the 
plateau suddenly appears above H c , it should be related 
to the specific scattering mechanism in the incommen- 
surate phase, i.e., the scattering by the domain- walls. 
Therefore, the field independence in the I-AF phase sug- 
gests that the soliton population and their distribution 
are not affected by the magnitude of the applied mag- 
netic field, indicating that the solitons are frozen. Since 
the plateau in n(H) is present only below Tjv [see Fig. 
3(b)], it is the specific feature of the I-AF phase, and the 
AF order must be essential for the soliton freezing. 

In the I-AF phase, it is suggested that solitons are 
preferably distributed around the impurity sites and that 
staggered-|moments are distributed also in between these 
solitons. Rp Since the spin-flop transition takes place at 
~ 1 T,E3Ey the staggered moments are directed almost 
perpendicular to H in the I-AF phase. When the mag- 
netic field increases, there are two ways to gain the mag- 
netic energy AMH, where AM is the additional mag- 
netization: one is to append other solitons in between 
the impurity sites, though this position is not so favored 
compared to the impurity sites, and the other is to cant 
the staggered moments to the direction of H . The sec- 
ond process can give an origin of the soliton "freezing" in- 
duced by the antiferromagnetism, because the soliton dis- 
tribution is preserved under further magnetic field. Al- 
though the above scenario gives a plausible explanation, 
independent measurements by neutron scattering and/or 
NMR are required to elucidate this microscopic mecha- 
nism of the antiferromagnetism-induced soliton freezing. 

The dotted line in Fig. 3(b) indicates that K is pro- 
portional to T 2 when the sample is cooled below Tjv. 
Since the specific heat have the temperature dependence 
of cx T 3 , the diffusivity should behave as oc 1/T, the typ- 
ical temperature dependence of K p h when phonons are 
scattered by dislocations of solids, like defects. The ob- 
servation shows that the domain walls of the spin solitons 
behave as rigid dislocations, which would be consistent 
with the above picture of the frozen solitons in the I-AF 
phase. 
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V. SUMMARY 



The results of the thermal conductivity measurement 
on both pure and Mg-doped CuGeC>3 , accompanied with 
the specific heat measurement for the same crystals, 
help to further understand the high-field incommensurate 
phase of the SP system, the spin-soliton phase. First, it 
turned out that the phonon transport in the I phase is 
governed by the domain-wall scattering due to the sig- 
nificant phonon-soliton interaction, and is strongly influ- 
enced by spatial distribution of the solitons. Since the 
kink features in H and T dependence caused by the ap- 
pearance (change in the distribution) of the solitons tell 
us the transition into (in) the I phase, the H — T phase 
diagram of the Mg-doped CuGeC>3 (where the CuO chain 
is directly disturbed), is established. Based on our result 
of k(H) above H c , it is suggested that the spin solitons 
are frozen in the I-AF phase with the help of the AF or- 
dering, providing a clue for the full understanding of the 
spin distribution in this unusual coexistence phase. 
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